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^ ■ ABSTRACT 

, Context. After leaving the main sequence, massive stars undergo complex evolution, which is still poorly understood. With a popula- 

tion of hundreds of OB stars, the starburst cluster Westerlund 1 offers an unparallelled environment to study their evolutionary tracks. 
Aims. We characterise a large sample of evolved OB stars in the cluster, with the aim of determining cluster parameters and place 
stars in an evolutionary sequence. 

, Methods. We used the F0RS2 instrument on the VCT to obtain intermediate-resolution spectroscopy over the range 5800-9000A of 

about a hundred stars selected as likely members of the cluster based on their photometry. We developed criteria for their spectral 
classification using only spectral features in the range observed. We discuss these criteria, useful for spectral classification of early- 
' type stars in the GAIA spectral region, in the appendix. Using these criteria, we obtain spectral classifications, probably accurate to 

■ one subtype, for 57 objects, most of which had no previous classification or a generic classification. 

Results. We identify more than 50 objects as OB supergiants. We find almost 30 luminous early-B supergiants and a number of less 
luminous late-0 supergiants. In addition, we find a few mid B supergiants with very high luminosity, some of them displaying signs of 
heavy mass loss. All these stars form a sequence compatible with theoretical evolutionary tracks. In addition, two early B supergiants 
also show indication of heavy mass loss and may represent the evolutionary phase immediately prior to the Wolf-Rayet stage. We 
Q . investigate cluster properties using the spectral types and existing photometry. We find that the reddening law to the cluster does not 

i—{ • deviate strongly from standard, even though extinction is quite variable, with an average value Ay = 10.8. Though evolutionary tracks 

' for high-mass stars are subject to large uncertainties, our data support an age of S 5 Myr and a distance » 5 kpc for Westerlund I. 

, Conclusions. The spectral types observed are compatible with a single burst of star formation (the age range is very unlikely to be 

> I Myr). Westerlund 1 shows its potentiality as a laboratory for massive star evolution, which can be fulfilled by detailed study of 
the population presented here. 

^ ' Key words, stars: evolution - stars: early-type - - open clusters and associations: individual: Westerlund I 

o 

1. Introduction HR diagram of massive stars during the H-shell burning phase 

. (Humphreys & Davidson 1994; van Genderen 2001). 

rn Massive stars play a crucial role in the dynamical and chem- ^here is general agreement that the WR stage corresponds to 

O , ical evolution of galaxies, providing a major source of ion- jj^.^^^g ^^^^i^g ^^^^^^ f^^^ ^^^^ ^j^^ ^^^y 

2 ■ ising UV radiation mechamcal energy, and chemical enrich- ^ot implies that massive stars do actually loop bluewards. But 

; ment (e.g., Massei|2003|; |van der Hucht et al. || 2003|) . As they ^j^^^^ ^ complex zoo of transitional objects, comprising blue 

> , evolve from the main sequence (MS) towards the Wolf-Rayet supergiants (BSGs), red supergiants (RSGs), yellow hypergiants 

(WR) phase, massive stars must shed most of their outer lay- (yhGs), LBVs and OBfpe/WNVL stars, whose identification 

^ ers. Models predict that massive stars will evolve redwards at ^^^^ particular evolutionary phase is a matter of guesswork. 

^ approximately constant Ltd, but there is an apparent dearth of Understanding this evolution is, however, crucial because the 

■ ■ ■ yellow supergiants with Lboi comparable to the brightest hot ^^^^ j^^g j^ring this phase completely determines the contri- 

stars. The lack of yellow hypergiants above a certain luminos- ^^^^^^ ^^at the star will make to the chemistry of the ISM and 

ity defines the observational Humphreys-Davidson (HD) limit, ^^^^ ^^^^ post-supernova remnant it will leave, 

which seems to imply that stars hit some kind of instability when , ^ ^ , . , , . , . 

they reach this area of the HR diagram, and then lose mass at a Unfortunately massive stars are scarce and, as this phase is 

formidable rate. This phase of enhanced mass loss is generally ^^'"3' .^^ort on evolutionary terms, examples of massive stars in 

identified with the luminous blue variable (LBV) stage, but there ti-ansition are rare. For most of them, distances are unknown and 

is no cleai- understanding of the actual evolution through the luminosities can be determined at best to order-of-magnitude 

accuracy, with the uncertainty feeding through to other parame- 
ters (Mt, Rt). As a consequence, the difficulty in placing these 

T , „ . , TZ , objects within an evolutionary sequence is obvious. 

Send ojfprmt requests to: I. Negueruela ■' ^ 

* Based on observations collected at the European Southern The starburst cluster Westerlund 1 (Wd 1) may represent 

Observatory, Paranal, Chile (ESO 71.D-0151, 073.D-0327, 081.D- a unique lab oratory for addressing these issues. With a mass 

0324) ~ 10^ Mq (iBrandner et al.l 120081: iMengel & Tacconi-GarmanI 
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I2OO8I) . it is young enough to contain a large number of OB stars, 
but old enough to contain stars in all the evolutionary stages. 
Large populations of WR stars and transitio nal massive stars 
(iNegueruela & ClarkI 120051: IClarket al.ll2005l henceforth COS) 
have been identified. In this paper, we set their evolutionary con- 
text by characterising the population of moderately evolved mas- 
sive stars in the cluster, i.e., objects that have started to move 
away from the MS, but have not yet reached the region of insta- 
bility. With this information, we are able to constrain the cluster 
parameters and check the agreement with evolutionary tracks, 
which supports the idea of an (approximately) single-age popu- 
lation. 



2. Observations 

Observations of stars in Wd 1 were carried out on the nights 
of 2004 June 12th and 13th using the spectro-imager FORS2 
on Unit 1 of the VLT (Antu) in three diff'erent modes: longslit, 
multi-object spectroscopy with masks (MXU) and multi-object 
spectroscopy with movable slitlets (MOS). We used grisms 
G1200R and G1028z, which provide the highest dispersions, 
to obtain intermediate -resolution spectroscopy. G1200R has a 
nominal dispersion of 0.38 A/pixel over the spectral range 5750- 
73 lOA. G1028z has a nominal dispersion of 0.42 A/pixel over the 
spectral range 7730-9480A. Together, they cover most of the red 
and nearest infrared range. 

For the longslit mode, we used a 0'.'3 slit to obtain the highest 
resolution available. Measurements of arc lines give FWHM of 
~0.9A with the G1200R, which results in R slightly above ~ 
7000. With the G1028z, we measure FWHM- 1.3A, implying 
a resolution slightly below ~ 7000. For the multi-object modes, 
we used 1" slits (10" long in the MOS mode and of variable - 
< 10" - length in the MXU mode). This results in a resolving 
power R = 2200 in the red range (the exact wavelength coverage 
depends on the position in the CCD) and R = 2500 in the near-IR 
range. 




Wavelength (A) 



Fig. 1. /-band spectra of two of the visually brightest B s uper- 
giants in Wd 1, tentatively classified as hypergiants in |C05L 
The top spectrum corresponds to W42a, which we classify as 
B9 Ia'^\ based on the features presented here. Note that the Ca 11 
triplet lines are already comparable in strength to Pa 13, 15 & 
16. The bottom spectrum corresponds to W33, B5 Ia*^\ Note the 
prominent P-Cygni profile in O i 8446A and the emission fea- 
tures at 8913 and 8927 A, seen in the spectra of all the mid-B 
supergiants in Wd 1 . 



We also used the MOS mode with grism G150I to obtain 
low-resolution spectroscopy of stars in the field. This mode gives 
a resolving power R - 260 over a range ~ 6000-1 lOOOA (the ex- 
act range depending again on the position within the CCD). This 
mode was used mainly for Wolf-Rayet stars, but other objects 
were also observed. 

The field of view of F0RS2 is 6.'8 x 6.'8 . Within the central 
5^0 X 5.'0 field of view, we selected our targets from the list of 
likely members of lC05l . For the external regions, targets were se- 
lected at random amongst relatively bright stars. In total, we took 
three MXU and one MOS mask with both G1200R and G1028z 
grisms, two further MXU masks with only the G1200R (these 
were aimed at relatively faint objects, which were expected to 
be OB stars near the MS and so not to have strong features in the 
range covered by G1028z) and three MOS masks with G150I. 
This resulted in ~ 120 stars observed with G1200R, ~ 75 stars 
observed with G1028z and ~ 50 stars observed with G150I. 
More than 90% of the spectra turned out to correspond to OB 
stars, and hence cluster members. 

In this paper, we restrict our analysis to the brightest OB 
members of the cluster, leaving analysis of the unevolved popu- 
lation of OB stars for a future work, where our FORS2 data will 
be combined with FLAMES and ISAAC observations. We com- 
plement our dataset with observations of a few additional bright 
OB stars which could not be observed during the FORS2 run. A 
few stai-s (W5, W7 & W52) were observed with the ESO Multi- 
Mode Instrument (EMMI) on the New Technology Telescope 
(NTT) at the La Silla Observatory (Chile) on the nights of June 
5-8 2003, using gratings #6 & #7, giving a coverage similar 
to that of th e FORS observations at comparable resolution (see 
IClarketalJ l2010 for details). Some other stars were observed 
with the NTT and EMMI on the nights of 15 - 17 February 2006. 
On this occasion, we used grism #6, which covers the range 
5800-8650A at a resolution R ~ 1500, though several objects 
are only detected from ^ 6500A. 

For completeness, we include a few other objects which have 
been observed only in the / band with the multi-object fibre spec- 
trograph FLAMES^_on_flie VIXDetails of those observations are 
reported in ' Ritchie et al.l (l2009ah . The spectra cover the range 
8484-9000A, with a resolving power R ^ \6, 000. 



3. Results 

We have selected for our sample those bright stars not 
listed as yellow or red supergiants (C05j) or Wolf-Rayet stars 
jCrowther et al.ll2006al) . The sample consists of 57 stars spectro- 
scopically identified as luminous OB cluster members. They are 
listed in Table 1 . 

Wd 1 is too obscured to allow the acquisition of classification 
spectra. Even with 8-m class telescopes, moderate SNR cannot 
be achieved at wavelengths shorter than 5500A with reasonable 
exposure times. We have been forced to rely on the data avail- 
able, covering the range between ~6000A and 9000A. Criteria 
for classification in this range were first explored in iC05L and are 
analysed and developed in Appendix A. 

As demonstrated in the appendix (see also Appendix A of 
C05), the classification of supergiants of spectral types later than 
B2 can be achieved with rather high accuracy, because of the 
presence of abundant metallic lines. For earlier spectral types, 
classification criteria are scarcer and most of them are sensitive 
to both spectral type and luminosity. In most cases, by using sev- 
eral of them at the same time, we can narrow the spectral range 
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Wavelength (A) Wavelength (A) Wavelength (A) 



Fig. 2. Ha spectra of the most luminous B-type stars in Wd 1 and som e com parison stars. The left panel shows the spectra of 
the three luminous B-type objects tentatively classified as hypergiants in IC05l The middle panel shows two objects with spectral 
types around BO with evidence for heavy mass loss. The right panel shows, as comparison, the spectra of the hypergiant candidate 
HD 80077, the hypergiant HD 190603 and the supergiant HD 148688, which has the strongest emission feature amongst all the 
objects examined in Appendix A. All the spectra have had the continuum normalised and are offset for display. 



for a given spectrum generally down to half a spectral subtype. 
This is achieved here by combining the red and /-band spectra. 

There are three very luminous B supergiants in Wd 1, which 
were tentatively classified as hypergiants in ,CQ5l The /-band 
spectra of W33 and W42a are displayed in Fig. [T] Based on 
them, they can be confidently classified as B5 la'^' and B9 la*^'. 
W7 is very similar to W33, and was also classified B5 la^^'Q. 
These objects display weak P-Cygni profiles in the Oi 8446A 
line, and much stronger P-Cygni profiles in Ha (see Fig. |2]i. 
Another interesting feature in the red spectra of W7 and W33 
is the C 11 6578, 6582A doublet. These lines are hardly seen in 
MS B-type stars, but appear weakly in B2-B5 stars of moderate 
luminosity and can be pro minent in supe rgiants. Its maximum 
strength occurs around B3 CWalbornlll980l) . 

Amongst Wdl stars, the strongest Cii is seen in t he spe c- 
tra of W70 and W71, both of which were classified in |C05| as 
B3 la based on their /-band features (strong O i 8446A line by 
the side of Pa 16, presence of Ni lines around Pa 11). Another 
object with prominent C ii doublet is W57a. From the ratio of 
Oi 8446 A to Pa 16 and the intensity of the Hei lines, this ob- 
ject is either later or more luminous than W70 (Fig. O. From 
the colours and magnitudes of these stars, the last possibility 
seems unlikely (all three have about the same magnitude in /, 
while W57a is less reddened). This is supported by the weaker 
Cn doublet. W57a shows distinctive Sin 6347 & 6371 A, as well 
as Ne I 6402A. These lines are also strong in W70, but weaker 
in W71. Therefore, we take B4Ia for W57a and B3 la for W70. 
W71 is slightly earlier, and we classify it B2.5 la. 



' Note that lCIark et alj l l201Cl) observe pulsational variability in these 
objects and speculate that they may induce variations of up to 1 spectral 
subtype 



Then we have W2a and W28, presenting moderately strong 
Cii 6578,6582A (Fig. gll and still a conspicuous Oi 7774A 
triplet, a strong feature in the later supergiants. The Oi 8446 A 
line is seen j ust as an inflection on the side of Pa 1 8 (Fig . O . Both 
display strong P-Cygni profiles in Ha and strong N ii 6482A. 
The presence of this line is c onsidered a signature of la class 
supergiants bv iLennonI (Il993h . We classify these objects B2Ia. 
Wll is very similar and receives the same classification. W23a 
has a very similar /-band spectrum. It also shows N n 6482A. 
However, Oi 7774A is missing, and the Cii doublet is very 
weak. These features are difficult to reconcile with any spectral 
type. As this star is very bright, we suspect that the spectrum 
may include an earlier type companion. W8b and W52 look 
slightly earlier than this previous group, but still have strong 
Cii 6578,6582A and perhaps Oi 8446A on the wing of Pa 18. 
We classify them as B 1.5 la. 

W78, still showing moderately strong C n doublet lines, has 
weaker O i 111 A A. The absorption trough in the P-Cygni profile 
of Ha has almost disappeared in this object. The spectrum of 
W78 extends further bluewards than most others, and we can 
identify a very strong Sim 5740A line, confirming a spectral 
type not far away from B 1 . Similar features are seen in W46a, 
W43b and W19. Though there is some evidence that these last 
three may be very slightly earlier, we classify all these stars as 
Blla. 

W18 and W21 show inconspicuous Cii absorption lines and 
broader Paschen lines, though we can still see them up to Pa 20. 
The He i lines in the /-band are now almost vanished except for 
the line at 8780A, still moderately strong. In W21, we know that 
it does not show conspicuous Si iii 5740, confirming that it is 
earlier than B 1 . Another similar spectrum is W6 1 a, showing pure 
emission in Ho'. C iii, if present, is very weak. We classify these 
three objects as BO. 5 la. 
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Fig. 3. Sequence of /-band spectra of blue supergiants in Wd 1, showing the evolution of the main features. Note the disappearance 
of N I features between Pa 12 and Pa 13 and 1 8446A around B2 (compare to Fig. [T] where lines are identified) and the development 
of Pa 16, which shows the presence of a strong C ra line for stars BO and earlier. 



W43a shows pure emission in Ho-. The Paschen lines are 
shallower and we can only see them up to Pa 18. Cm 8500A 
is starting to contribute strongly to Pa 16. There is no sign of 
Cn or Nil absorption now, but strong Cn wind emission lines. 
This suggests a spectral type BO la. This object is a single-lined 
spectroscopic binary (Ritchie et al. 2009a). W34 has a similar 
spectrum and receives the same classification. 

W232 and W54 have broader Paschen lines, even though we 
see up to Pa 18. W54 presents a neutralised Ha, while W232 has 
it weakly in absorption. It also shows a very weak Si m 5740A 
line. The lack of Ha emission likely indicates a lower mass 
loss and hence luminosity. We classify these objects as BO. 5 lab. 
Another similar spectr um is W6a, which shows variable Ha and 
hard X-ray emission (IClark et al.ll2008h . and is a likely binary 



(see below). As discussed in Appendix A, the strength of Ha 
emission correlates well with luminosity, but there are many ex- 
ceptions. 

W60 and W63a have similar spectra, with Ha neutralised or 
weakly in absorption, but show clear evidence of C iii 8500A. 
Paschen lines are broad, and we classify these two objects as 
BO lab. Based on the overall similarity, we classify W49, for 
which we only have a FLAMES /-band spectrum, as BO lab, as 
well. 

W61b also has Ha weakly in absorption, but has a prominent 
Cm line, while Paschen lines are visible up to Pa 18. We clas- 
sify this object as 09.5 lab, though a higher luminosity cannot 
be ruled out. Two other stars look like high-luminosity O-type 
supergiants. W74 has a very prominent Cm 8500A line, and 
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6600 
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Fig. 4. Sequence of red spectra of blue supergiants in Wd 1 . The 
shape of Ha changes from a P-Cygni profile to pure emission 
around spectral type BO. 5. The N n 6482A line is only prominent 
in la supergiants and its strength peaks sharply at B2. 



Paschen lines up to Pa 18. Its Ha line is weakly in absorption. 
We classify it as 09.5 lab, though it also could be slightly earlier 
or slightly more luminous. W24 is very similar, with Ha more 
clearly in absorption, and its Paschen lines show that it is either 
less luminous or earlier t han the other two. This object shows 
some spectral variability dClark et al.ll2010h and needs further 
study. We provisionally classify it as 09 1. 

Similar spectra, with Ho- close to being neutralised, are 
shown by W35 and W41 . As they have stronger C m 8500A and 
Pa 18 is not clearly detected, we classify them as 09 lab. Based 
on overall similarity, we also classify as 09 lab a number of ob- 
jects with lower-resolution (NTT) or incomplete (only /-band) 
spectra, W17, W25 and W38, though these classifications are 
less secure. W27 has a similar spectrum, but shows several indi- 
cations of being a binary (complex Ha, shallow broad Paschen 
lines, likely double He i 6678 A), in accord with its high X-ray 
luminosity dClark et al.ll2008h . 

W30a displays strong variability and clear signs of binarity 
(see also Clark etal. 2008). Some of the spectra suggest that 
Hen 8236A could be present, while Cm 8500A changes be- 
tween spectra. The Paschen lines are broad and shallow, typi- 
cal of a spectroscopic binary. This system contains two O-type 
supergiants, and one of them could be earlier than 09. 

For W3 1 , we only have spectroscopy at one epoch, but the 
double-peaked structure of Ha strongly suggests that this BO I 
supergiant has a companion, rendering the exact luminosity clas- 
sification impossible. WIO, as shown in Clark et al. (2008), is a 
double-lined spectroscopic binary, and therefore must contain 
two supergiants. 

To the South of the cluster, W238 does not seem to fit ex- 
actly with the spectral sequence. It has C n in absorption and Ha 
partially neutralised, suggesting a spectral type Bl lab or Bl lb. 
However, /18502A looks stronger than expected, and this could 
also be a composite spectrum. . 



All the other objects have Ha in absorption and a smaller 
number of Paschen lines. Based on the relative strength of the 
C III 8500A feature, we classify them 09-09.5 lb. We consider 
that having Ha in absorption to the point that it is stronger than 
the He i 6678 A line is the defining characteristic of luminosity 
class lb. However, we caution that sky subtraction is difficult in 
the very crowded regions close to the cluster core, and Ha could 
be affected by the sky. 

Taking all these factors into account, we estimate that the 
spectral types we derive are likely accurate to within half a sub- 
type for B supergiants. Luminosity classes are also likely to 
be correct to within one subclass. Types for O stars are not so 
well determined, though we expect to be correct within one sub- 
type. Luminosity classes are less reliable, as we lack proper cri- 
teria. We have to resort to line widths and shapes, which are 
also dependent on rotational velocity. As the rotational veloci- 
ties of O supergiants may be relatively high, there is a degree 
of degeneracy between luminosity and rotation. Also, the types 
derived depend strongly on the ratio of the C m 8500A line to 
neighbouring Paschen lines. This indicator is obviously sensi- 
tive to CNO anomalies. However, we do not find any spectrum 
where the strength of this line is openly in contradiction with any 
other criteria. Based on this, we believe that any CNO anomalies 
present in our stars are, at most, moderate. In the case of moder- 
ate C enhancement, a star classified as 09 lab might actually be 
a BO lb with enhanced C m. In the case of moderate N enhance- 
ment (and, hence, C depletion), a star classified BO-0.5 lab or lb 
could be an 09 la with a very weak C iii. None of these possible 
misclassifications looks too severe. 

Finally, Table 1 lists the spectral types of 55 OB supergiants, 
and two less luminous objects, in Wd 1 

3.1. Bright B supergiants 

Two B5 la and one B9 la supergiants show characteristics that 
lead us to classify them as hypergiants. For a start, they are 
> 1 mag brighter (after correcting for extinction; see below) 
than any other B-type supergiant. In addition, they show signs of 
heavy mass loss, such as strong P-Cygni profiles in Ha (Fig.|2l) 
and the presence of a P-Cygni profile in Oi 8446A (Fig. [TJ. 
IClark et al.l (1201 Ol) present evidence for variability in these three 
objects. At this stage, we do not intend their classification as 
hypergiants to have any implication about their evolutionary 
status. Most of the B5-9I supergiants that have been subject 
to detailed analysis fall on tracks corresponding to stars with 
< 25 Mq jMcErlean et al.l 1 1999; Markova&Puls 2008), and 
hence it seems natural to find that the late B supergiants in Wd 1, 
most likely descendants of more massive progenitors, are more 
luminous. Analysis of their chemical compositions will cast light 
on their evolutionary status. 



3.2. Transitional objects with emission lines 

Two objects displaying strong Ha emission and Hei emission 
can be classified as early B objects, W5 and W13. Their Ha 
pro files can be see n in Fig.|2](wider range spectra are displayed 
in lClark et al.l2 008). Their /-band spectra are displayed in Fig.|5] 
where they are compared to the WN9 star W44 (= WR L). 

W5 displays very intense Ha emission, P-Cygni profiles in 
He I AA 5785, 6678 & 7065, C n wind features (not only the 723 1, 
7236A doublet is in emission, but also 6578, 6583A) and some 
N II lines weakly in emission. In contrast, its / band spectrum 
lacks any emission features and seems typical of a supergiant 
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Wavcleneth (A) 



Fig. 5. /-band spectra of three transitional objects with emission 
Hnes. W44 (=WR L) is a late Wolf-Rayet star, presenting emis- 
sion components in all Paschen lines. W13 is an eclipsing binary. 
One of the components is the emission-line object. W5 shows no 
emission lines in the /-band in spite of strong emission features 
in the red (see Fig.|2] 

around B0.5 (though weak emission components in the Paschen 
lines may not be seen at our resolution). In view of this, we clas- 
sify this object as a BO. 5 la^ star, where the hypergiant classifica- 
tion is due to the several indicators of heavy mass loss. However, 
we note that this object is intrinsically less luminous than other 
stars classified as BO. 5 la. In view of this and its similarities to 
WN9, it is tempting to identify W5 with an evolutionary phase 
immediately prior to the WR stage. Again, a careful analysis of 
its chemical composition may provide further clues as to its evo- 
lutionary status. 

W13 has a similar spectrum, but analysis is at present com- 
plicated by the fact that t his object is a double-lined eclips- 
ing binary ( Bonanosll2007l: [Ritchie el al. 2009a). Pending a de- 
tailed analysis of the system, the data seems consistent with an 
emission-line component similar to W5 and an absorption-line 
object. Though the Pa 16 featu re in Fig. 5 is very prominent, 
several FLAMES spectra (cf. Ri tchie et al. ''2009a) show moder- 
ately strong He i lines, indicating that the companion is likely an 
early-B supergiant. 

3.3. Likely binaries 

Several spectra show evidence for binarity. In WIO we see dou- 
ble He I lines. In W23a, we see features incompatible with a sin- 
gle spectral type and variations between two epochs. W30a and 
W31 show evidence of mass transfer in the form of complex 
and (in W30a) variable Ha. W30a and W27 show shallow lines, 
typical of double-lined bina ries. W6a also shows important vari- 
ations in the shape of Hor dClark et a l. 2008) and its FLAMES 
spectra show very shallow broad lines. The spectra of W27 have 
lower resolution than those of other sources, but the shape of 
the He i lines is also suggestive of two components, while Ha 
shows a complex emission profile. Of these, W6a, WIO, W27 
and W30a are hard X-ray sources (Clark et al. 2008), the latter 
being the most luminous and hardest det ection in the cluster after 
the magnetar CXO J164710.2-455216 (IClark et al.ir2008.) . W27 
is also a relatively bright X-ray source (Lx ~ 10^^ erg s"^). 

Considering the relative brightnesses of OB supergiants and 
MS stars, every object showing double lines must be considered 
to contain two supergiants. Given the relative shortness of this 
evolutionary phase, the presence of at least two (and likely more) 




8500 8600 8700 8800 

Wavelength (A) 



Fig. 6. Sequence of FLAMES/GIRAFFE /-band spectra of some 
of the earlier spectral types observed in our sample. The bottom 
two spectra have lower SNR, but the Paschen lines are clearly 
seen to have much broader wings than the more luminous ob- 
jects and display much weaker C iii 8500A. Note the difficulty 
in defining the continuum around the area of the DIBs at 8620A 
and 8648A, in spite of the high resolution. The spikes seen in 
some spectra are badly removed sky Unes. 

objects displaying the spectra of two supergiants suggests that 
there is a high number of systems which actually are similar- 
mass binaries, in good agreement with the preliminary results of 
a radial velocity surve y, designed to determine the binary frac- 
tion in Wd 1 (.Ritchie et'al.ll2009ah . 

3.4. The less luminous objects 

In order to illustrate the characteristics of O-type stars in the /- 
band, we have included in our sample two stars with significantly 
fainter / magnitude. These are W6b and W50b. Their spectra are 
displayed in Fig.|6l compared to some other objects that were ob- 
served only with FLAMES. The small spectral coverage makes 
all classifications somewhat uncertain, but Fig. |6] shows that a 
few spectral features may still provide some guidance. The main 
discriminant for luminosity class is the shape and depth of the 
Paschen lines. 

The two faint stars have significantly broader Paschen lines 
than any star that we classify as supergiant. The limit between 
the wings of the lines and the continuum is difficult to define, 
in opposition to the sharply defined lines of the supergiants. 
Lacking any further criteria, we tentatively classify these two 
stars as luminosity class III, noting that they may well be of 
lower luminosity - as suggested by their dereddened magnitudes 
(see below). They are representative of the large population of 
unevolved massive stars present in Wd 1, which we will study, 
with a more extensive dataset, in a future paper. 

4. Discussion 

We have spectroscopically identified 55 OB supergiants in Wd 1 . 
The sample is still very far from complete. On the one hand, the 
extended halo around the cluster is not fully explored. On the 
other hand, many stars in the cluster core, with magnitudes and 
colours similar to the O-type supergiants Usted here, still lack 
spectroscopy. 
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Fig. 7. Correlation between the E{V - I) and E{R - I) colour 
excesses for OB supergiants in Wd 1 . Individual excesses have 
been derived from the photometry listed in Table 1 and the in- 
trinsic colour calibration. 



The significance of this huge population has to be discussed 
within the context of their membership in a single cluster, Wd 1 , 
which hosts a large population of massive stars in more advanced 
evolutionary stages. Th e cluster co ntains one probable active hot 
LBV (W9; Doughertv et al.ll2010h. 7 A/F hy pergiants (includ- 
ing the cool LBV W243nRitchie et al.ll2009bh . 4 very luminous 
RSGs ( Clark et al. 201 Oh and a large population of WR stars 
("Crowther et al. 2006a). The OB supergiants provide us with the 
possibility to constrain the evolutionary paths of these evolved 
stars. Unfortunately, the parameters of Wd 1 are still poorly de- 
termined, even if significant progress has been made in recent 
years. The data presented here can set strong constraints on many 
of these parameters. 



4.1. Extinction 

The distance to the cluster is, at present, not very well con- 
strained. Photometry is affected by very strong reddening. 
Analysis of the E{B - V) colours suggested that the reddening 
could deviate from the standard law (COS), with an extinction 
Ay ~ 12 mag. With the new accurate spectral types, we can con- 
duct a deeper investigation of the reddening to the cluster 

We have taken all the objects with photometry in Table 1 
and calculated their colour excesses, making use of the intrin- 
sic colours for OB supergiants (Wegner 1994). Figure |7] shows 
the run of E(R - I) against E(V - I). The two colours cover a 
wide range, showing that the reddening is very variable across 
the face of the cluster E{V - I) goes from ~ 5.0 mag to 
~ 6.4 mag, with a single object, W84, showing rather lower red- 
dening E{V - /) = 4.5. The two colours are very well correlated 
{R = 0.99 for a fit without independent term). The best fit gives 
E(V - I)/E{R - /) = 1.97 + 0.01, which c ompares very well 
with t he standard reddening law (e.g., 1.95 in lRieke & Lebofskvl 
11985 V The average value is E(V - /) = 5.5 + 0.3 (1-cr standard 
deviation). Similarly good fits are obtained for E{V-I)/E(V-R) 
and E(R - I)/E(V - R) (which, obviously, are not independent), 
giving values very close to the standard reddening law. 

This suggests that the extinction is close to standard, in spite 
of the difficulty with the E{B - V) colours. The correlation be- 
tween E{B - V) and E(V - I) is poorer If we force the fit to go 



through the origin, we obtain E(V - I)/E(B - V) = 1.31 + 0.01, 
with R = 0.97. This value is almost identical to the average 
of the ratios E(V - I)IE{B - V) for all the stars, 1.34 + 0.15. 
The value is not compati ble with the standard value, (e.g., 1.6 in 
iRieke & Lebofekvl[T98l . but, because of the large dispersion in 
individual values, is less than 2cr away. In view of this, we sus- 
pect that the difficulty stems from unreliable values for_B, per- 
haps a zero-point offset. Even though the photometry in|C05|was 
very carefully calibrated, stars in Wd 1 are much redder than any 
standards used, and have very faint B magnitudes. In any event, 
the average value directly determined from the observed {B - V) 
colours is E{B -V)- 4.2 ± 0.4 (again, 1-cr standard deviation). 

Another reliable indicator of interstellar extinction is the 
presence of Diffuse Interstellar Bands (DIBs). Unfortunately, 
the correlations between the Equivalent Width (EW) of most 
bands and reddening saturates at relativ ely low value s (e.g. , 
|Munai-i & Zwitter 19971; ICox et alJl2005l) . iMunari et alJ (l2008h 
find that the 8620A DIB is an excellent tracer of reddening, with 
a very tight coiTelation E{B - V) = 2.72(+0.03) x EW (A). 
As discussed in the appendix, measurement of the intensity of 
this band is not easy at our resolution. However, we may expect 
the high number of objects to compensate for the uncertainty 
in individual values. Unfortunately, there is no obvious correla- 
tion between EW(DIB) and E{V - I) for our stars. Objects in 
Wd 1 show an average value EW(DIB)= l.OOA, with a disper- 
sion cr - 0.09 and individual values ranging from 0.85 to 1.15, 
with no correlation to reddening or spatial position. 

If we restrict the search to objects with high-resolution 
FLAMES spectra, we find a similar result. The average is now 
1.1 1 A, showing that the lower resolution data tend to give lower 
values because of the poorer definition of the continuum. The 
dispersion is only cr - 0.08, but still no correlation with E{V-I) 
or location is obvious. 

This suggests that this DIB is saturated, at least along this 
line of sight, at the reddening of Wd 1, and the dispersion 
amongst Wd 1 stars is due simply to the difficulty in determining 
the continuum. As discussed in the appendix, this area contains 
a broad shallow DIB centred on 8648 A (cf. lMunari et al.ll2008l) . 
which becomes compo unded wi th some weak absorption lines. 
The correlation of Mu nari et al.l (2008) would give an average 
E(B - y) = 3.00, which is smaller than both the observed value 
and a direct extrapolation from E(V - R) and E(V - 1) under the 
assumption of a standard reddening law. 

There is another sharp DIB in our range, that at 6614A. 
Unfortunately, its corr elation to E(B - V) seems to be poor at 
high reddening values ( Cox et alJ2005h . Again, we find no obvi- 
ous con-elation to E(V - I). Most stars show EW(DIB)= 0.60A, 
with standard deviation cr = 0.05. As this value i s comp arable 
to the highest EW(DIB) measured by ICox et al.1 (|2005') along 
sightlines with high extinction, it is possible that this line is also 
saturated, though measurements at higher resolution should be 
desirable to check this hypothesis. 

The only DIB that has been calibrated against E{B - V) and 
might perhaps not be saturated at the reddening of Wd 1 is the 
narrow band at 5797A. Unfortunately, this is just outside the 
range of our spectra. It is seen in a few of the MXU spectra, 
due to their variable spectral range, but the continuum cannot be 
defined well. 



4.2. Distribution of spectral types 

The earliest supergiants in Wd 1 seem to have spectral types 
around 09. This is defined by the presence of strong C ni 8500A 
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(while the neighbouring Paschen hnes are weak or absent) and 
lack of Hen 8238A. Even though our spectra have not been 
corrected for telluric absorption, humidity was very low on the 
night of June 13th, and the stronger Hen line of earlier spectral 
types should be noticeable. None of the stars classified as 09 
or 09.5 seems sufficiently luminous to receive a la luminosity 
class, based on their relatively broad Paschen lines and lack of 
detection beyond Pa 17 (but see below). W74 (where Pa 18 is 
weakly seen), is the most luminous 0-type supergiant based on 
these criteria. 

There are a fair number of luminous (la) supergiants cover- 
ing the range B0.5-B2. Later spectral types become less numer- 
ous, with one each of B2.5 la, B3 la and B4 la, and two B5 la*^^ 
bright supergiants. After this, we only have a B9 la'^' object be- 
fore we reach the A hypergiants (w hich include the present phase 
of the LBV W243: iRitchie et al.ll200 9b). This distribution hints 
at some sort of break around spectral type B2-3, perhaps sug- 
gesting that stars in this mass range finish their H-core burning 
phase as ~B2 supergiants. 

Though our sample of OB supergiants is far from complete, 
there have been no systematics in the choice of targets. The 
level of incompleteness is difficult to assess, as photometry is 
not available for all the stars in the crowded cluster core area. 
We estimate that perhaps ~ 10 stars which appear bright enough 
in / to be supergiants have not been observed in the cluster core 
(most notably W12b, W39a & W40a), with a similar number 
in the outskirts (e.g., W53 or W56a). Of course, a much larger 
population of less luminous OB star remain to be explored. 

4.3. Mass 

The mass of Wd 1 has been estim ated in two different ways. 
iMengel & Tacconi-Garman 1 120081) have used the radial veloc- 
ity dispersion of the ten stars brightest in the infrared (cr - 
8.4 km s"') to estimate a mass of ~ 1.3 x 10^ Mq. This estimate 
needs to be confirmed, as it is subject to a large uncertainty. The 
YHGs s how large variations in their radial velo cities due to pul- 
sations ( iRitchie et al.ll2009atlciark et al.ll2010l) . and some of the 
stars used present important spectral pecuharities (e.g., W26 is 
associated with extended nebular emission). Because of this, and 
the assumption o f virial equilibrium, th is value should be taken 
as an upper limit. iBrandner et all (l2008l) . on the other hand, have 
used star counts in the infrared to set a lower limit on the cluster 
mass Mci ^ 5 x lO'* Mq. Again, there are important uncertainties 
coming into this estimation, such as the completeness correction. 

In particul ar, the pre-main-sequence isochrones used by 
IBrandner et al.l (|2008) indicate a significantly younger (t - 
3.3 + 0.2 Myr) and less distant (<i = 3.6 + 0.2 kpc) clust er. This 
results in very different masses for a given star. Recentlv. iNavlorl 
(|2009|) has shown that similar discrepancies are observed for 
a large number of young open clusters and suggested that this 
effect reflects a real deficiency in the calculation of pre-main- 
sequence ages. 

In any case, the integrated initial mass function, which is 
independent of the distance and age estimate, appears consis- 
tent with Salpeter's down to solar mass (IBrandner et al.ll2008h . 
The spectroscopic confirmat ion o f a large number of photomet- 
ric candidate members from ICOSi here indicates that the number 
of stars more massive than the main sequence turn-off is likely to 
be ~ 200. If we accept that the turn-off stars have M, ^ 25 Mq, 
a Kroupa IMF would then imply a mass ~ 10^ Mq, consistent 
with both determinations. 
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Spectral type 

Fig. 8. Observational HR diagram for Wd 1 . The absolute mag- 
nitudes. My, are estimated from the observed /, dereddened us- 
ing E{V - I) under the assumption of a standard law and cor- 
rected with th e tabulated (V - /)o for the spectral type from 
We gned (11994 . Circles represent objects classified as la super- 
giants; squares are lab classifications, diamonds are lb objects 
and triangles are stars with lower luminosity. 

An important source of uncertainty in both measurements 
is the role of binarity, which affects both radial velocities and 
transformation of magnitudes into masses. Different lines of ev- 
idence support a very high fraction of binaries with similar mass 
components amongst cluster members. Most WR stars show 
indirect indications of binarity: presence of dust in WC stars 
(Crowther et al. 2006a) and hard X-ray emission in WN stars 
(Cl ark et al.l 12008). Some of them also show photometric vari- 
ability and at leas t one is an ec l ipsing binary (.Bonanos,, 20071) . 
The pilot study bv IRitchie et al.l (l2009ah shows that the high bi- 
nary fraction extends to lower masses. 

4.4. Distance and age 

A recent determination of the distance to Wd 1, making use 
of atomic hydrogen in the direction to the cluster, gives d = 
3.9 + 0.7 kpc (Kothes & Dougherty 2007) , compatible with, 
though slight ly shorter than, estima tes based on the stellar pop- 
ulation (e.g.. Crowther et aD l2006a'). If we use the standard red- 
dening law to deredden the magnitudes and assume DM = 13.0, 
we find that luminosity class la stars typically have My between 
-6.4 and -6.90. This is not in bad agreement with existing cali- 
brations (e.g., Humphreys & McEkoy 1984), though the slightly 
longer distance of Crowther et al. (2006a) (DM ^ 13.4) would 
give a better agreement. The three mid-B supergiants W57a, 
W70 and W71 are definitively brighter, with -8 < My < -7. 
W7 and W33 are even brighter, with My -8.5 (which would 
become almost -9 for DM = 13.4), justifying their tentative 
classification as hypergiants. 

In Fig. [8j we plot an estimate of the intrinsic luminos- 
ity for all the objects with good photometry in Table 1 dis- 
tributed by spectral type and luminosity class. For this, we use 
E(V - I), under the ass umption of a standard reddening law 
(iRieke & Lebofskvl [19850 . to calculate A/, and then derive M; 
for a nominal DM - 13.5. Due to the uncertainty in E{B - V), 

^ Though a few are much fainter. In particular, W55 has -5.6, intro- 
ducing some doubts about its classification, wliich is based only on a 
FLAMES spectrum. 
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we prefer to perform our calculations with the red magnitudes. 
However, for ease of comparison, we transform to My simply 
by adding the intrinsic {V - I)o tabulated for the spectral type. 
The plot shows, in general, a good correlation between abso- 
lute magnitude and luminosity class, giving strong support to our 
classifications. The obvious exception are the objects classified 
as 09 SGs. These objects cover a broad span (> 2 mag) in My, 
present an overlap between 09Iab and 09Ib objects (though 
objects classified 09 lab are in most cases brighter than objects 
classified 09 lb) and have My's rather higher than expected (the 
09 lab objects are as bright as the BO la stars). 

This discrepancy may arise from the difficulty in separat- 
ing the effects of luminosity and fast rotation in O-type stars 
at this moderately low resolution. Some of the objects classi- 
fied as 09Iab may be fast-rotating 09Ia stars. Though plausi- 
ble, this interpretation presents two difficulties: first, all 09 la 
stars present should be fast rotators (perhaps not a serious objec- 
tion when we are talking of 3-4- stars); second, it would mean 
that we have a very large range of luminosity classes at spec- 
tral type 09 (from 09 III to 09 la), as opposed to other spec- 
tral subtypes. Alternatively, as discussed above, we can think 
that some of the objects classified as 09 lab are really later-type 
(~B0) stars with enhanced C iii 8500A, but this does not neces- 
sarily imply brighter My- Perhaps the simplest interpretation is 
that the objects we have classified as 09Iab are really 09.5Ia. 
As seen in Fig. 8, this small change in their classification would 
make their luminosities compatible with all the other members^. 

With this caveat, we may try to constrain the age of the clus- 
ter from the observed population. For this, we are forced to rely 
on the predictions of evolutionary models. As discussed above, 
the tracks followed by high-mass stars in the theoretical HR di- 
agram are subject to large uncertainties. As the exact mecha- 
nism for heavy mass loss is unknown, mass loss rates are nor- 
mally introduced by hand, following some recipe. Differences 
in mass loss rates, the treatment of rotation and the degree of 
convective overshooting assume d lead to rather different evo- 
lutionary tracks (e.g ., [Salasnich et alJll999t iMevnet & Maedej 
l2000t lEldridge et al]|20 08). Furthermore, the inclusion of bina- 
rity can substantially alter the evolution of massive stars (e.g., 
IVanbeveren etani2007l: lEldridge et alJl2008l) 

As discussed by ' Crowther et alj (l2006al) . single evolution- 
ary models have trouble reproducing the obs erved population of 
Wd 1. Geneva isochrones without rotation (iMevnet & Maede^ 
I2OOOI) predict that RSGs will be formed only for ages ^ 6 Myr, 
while even higher ages are required to produce RSGs with the 
high-rotation isochrones. On the other hand, the WR populations 
decline rapidly after 5 Myr The simultaneous presence of large 
numbers of cool super/hypergiants and WR stars is reproduced 
better by models including binaries, which also provide a better 
description of the ratio between WC and WN stars. Comparison 
of the number ratio between WR stars and cool hypergiants 
with iso chrones from Eldridge et al. (2008) led Crowtheret al] 
(l2006ah to propose an age of 4.5 or 5.0 Myr, with the progen- 
itors of the WR stars having initial masses in the 40-55 Mq 
range. 

Such age seems consistent with the observed distribution of 
blue supergiants. We observe O9-B0 supergiants with different 
luminosity subclasses, while at later spectral types we only see 
la supergiants. At ~ 4.5 Myr, these la supergiants should be de- 
scended from stars with initial masses Mini ~ 35 Mq. This is in 



^ As discussed in Section 3, many of the 09 lab classifications are 
less reliable than the rest, due to lower resolution or smaller spectral 
range observed. 



good agreement with pl acement on theoretical tracks of B la su- 
pergiants. For instance. ICrowther et al.l (l2006bh find that a sam- 
ple of BO-3 la Galactic supergiants lie between the tracks for 25 
and 40 Mq. 

As single star isochrones have been frequently used to date 
open clusters, here we will compare the properties of OB super- 
giants to Geneva isochrones. As the stars are still relatively un- 
evolved, the differences with respect to other single star tracks 
are minor, allowing easy comparison to the ages of other open 
clusters. 

In Fig. |9] we plot the same estimates of My derived above 
against Teft for all the la and lab supergiants. The effective 
temperatures u sed are an in t erpola tion between the observa- 
tional scale of iMartins et alj (l2005l) for late O -type stars and 
the temperature scale of ICrowther et al.l (l2006b h for B-type su- 
pergiants. We also plotjeveral isochrones from the models of 
iMevnet & Maedeil (1200 0). Taking into account that the My de- 
termination is subject to important uncertainties (for example, 
simply the rounding up of the second decimal in the determina- 
tion of the ratio between Aj and the colours amounts to a dif- 
ference > 0.1 mag) and the approximate character of the T^s 
scale, the general agreement is rather good. We note the follow- 
ing points: 

- The logf = 6.7 (5 Myr) isochrone without rotations would 
provide a good fit for a distance approaching 6 kpc. 

- The log t = 6.8 (6.3 Myr) isochrones (with and without rota- 
tion) provide a very good fit to the position of the B super- 
giants for distances ~ 4.5-5 kpc. These isochrones would, 
however, fail to reproduce the position of the O supergiants 
and almost certainly fail to reproduce the turn-off of the main 
sequence, unless all the objects we have classified as O-type 
supergiants are blue stragglers. 

- The location of the mid-B supergiants is compatible with 
their being evolvin g away fro m the main sequence. 

- The distance from ,Kothes & D oughertv (2007) (DM = 13.0) 
implies an age > 7 Myr for the cluster, which seems incom- 
patible with the population of WR stars. Their error bars, 
however, reach within the area favoured by our data. 

If we accept that a significant fraction of our stars are bi- 
naries containing two massive stars, they would be intrinsically 
brighter, allowing a longer distance (by ~ 0.5 mag) for a given 
age. 

Considering the uncertainties involved here, a definite an- 
swer will have to wait for accurate spectral classification of 
objects near the turn-off of the main sequence, which can be 
achieved with high quality infrared spectra. This will allow the 
anchoring of the infrared photometric sequence, and a much 
more accurate isochrone fit. The richness of the OB population 
in Wd 1 may then allow the creation of a more accurate tempera- 
ture scale based on fits of tailored synthetic spectra to high qual- 
ity spectra of our sample. Meanwhile, our data strongly favour 
d ^ 5 kpc and r k. 5 Myr for Wd 1 , an d seem to definitely 
rule out the values o f iBrandneret 10(120081) . based on pre-main- 
sequence tracks. As this paper contains a careful analysis of the 
data, we suspect that this r esult reflects in accuracies in the pre- 
MS isochrones, as found bv iNavioi1(l2009h for other young open 
clusters. 

4.5. The yellow hypergiants 

Wd 1 contains six luminous hypergiants wi th spectral types be- 
tween A5Ia+ and F8Ia+ (,Clark et al.ll201Q) . They represent a 
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Fig. 9. Semi-empirical HR diagram for bright OB supergiants 
in Wd 1 . The solid lines represent the Geneva isochrones with- 
out rotation for logf - 6.7 (5 Myr; top, blue) and logf - 6.8 
(6.3 Myr; bottom, brown). The dash-dotted line is the log t - 6.8 
isochrone for high initial rotational velocity. 



sizable fraction of the known Galactic YHGs. The evolution- 
ary status of YHGs is unclear un der current models. In several 
schemes (e.g.. ISmith et al.l [20041) . they are assumed to be stars 
that have already been RSGs and are coming back towards the 
blue in their way to become Wolf-Rayet stars. With the current 
dataset, the evolutionary status of the YHGs in Wd 1 is diffi- 
cult to assess. A simple estimation of their luminosities using 
the same procedure utilised for the OB supergiants suggests that 
their bolometric luminosities are not very different from those 
of the B hypergiants (and, in the case of W8a, substantially 
lower) Q. 

We have to be careful with this estimation, though. 
Spectroscopic monitoring of the YHG W265 has shown it to 
present strong pulsational variability, which leads to changes 
in spectr al type between F l la^ and F5 la^, with a period of 
~ 100 d (IClarketal .1120101) . There is a high likelihood that the 
other YHGs experiment similar variations. Such variability in 
spectral type should be reflected in significant changes in the in- 
trinsic colours and bolometric corrections. Our photometric and 
spectroscopic data are not simultaneous, and so the two datasets 
may correspond to different spectral types. Moreover, and per- 
haps more importantly, strong pulsations may lead to heavy mass 
loss, and the assumption of typical colours for A-F supergiants 
may not be appropriate if the stars suffer from intrinsic absorp- 
tion. 

Therefore a detailed spectroscopic analysis, including model 
fits for abundance determination, seems instrumental in deriv- 
ing the evolutionary status of the YHGs. So far, such an analy- 
sis has only been conducted for the LBV/A3 la^ W243, which 
shows very clear sig ns of heavy element pro cessing in its atmo- 
sphere. As noted by ICrowther et al.l (l2006al) . single star evolu- 
tionary models do not reproduce the simultaneous presence of 
luminous cool supergiants and WR stars in Wd 1. Models in- 
cludin g binary evol ution fare better in this respect but, as pointed 
out bv IClark et al] (l20To'). no current model seems able to ac- 
commodate the large population of cool supergiants at the es- 
timated age. Understanding the evolutionary relations between 



Note that we lack photometry for such interesting objects as W42a 
(B9Ia+), W32 (F5 Ia+) and W243 (LBV/A3 Ia+). 



the evolved stars in Wd 1 will therefore provide very strong con- 
straints on evolutionary models. 

4.6. Spatial distribution 

Based on star counts in the infrared. ferandner et al.l(l2008h found 
the cluster to be clearly elongated. The diff use X-ray emissio n 
from the cluster shows a similar elongation (iMuno et alJl2006l) . 
with the major axis at Position Angle 13° + 3°. 

In Fig. [To] we compare the distribution of the OB SG to that 
of the more evolved stars. For this latter set (comprising A/F hy- 
pergiants, M supergiants and WR stars), we can be certain of 
having a complete sample. The distributions of the two sets are 
not significantly different, revealing a very compact core, which 
contains more than half the objects of interest. This core mea- 
sures only ~ 60" x 45", with the centre approximately defined 
by the positions of W27, W29 and W30, and is elongated along 
the approximate NE-SW line. In addition, there is a secondary 
clump to the South, with a number of objects likely defining a 
more symmetrical halo. We note, however, that the division in 
two clumps is not obvious in eithe r the diffuse X- ray emission 
or the infrared star counts. Indeed. [Brandner et all ([2008) place 
the centre of the cluster in a position very close to WR L, in the 
area of lower density of supergiants. This may suggest that the 
division in two clumps is not real, but due to a random decrease 
in the number of optically bright stars. In such case, the shape 
of the cluster, as defined by its supergiant (and WR) populations 
would resemble more closely an arc than an elongated structure. 

Some objects, like W228b and W373, are located more than 
150" away from the putative core centre, hinting at the possibil- 
ity that Wd 1 is surrounded by a very extended halo. Though this 
possibility has not been explored yet, it is given strong support 
by the presence of 3 WR stars likely associated with the cluster 
at m uch higher distances (WR stars N, X and T in the notation 
of iNegueruel a & Clarkll2005HCrowther et al. 2006a). 

At this level of analysis, we find no strong indication of any 
difference in the relative distributions of evolved and unevolved 
stars, which could point to an age difference between regions in- 
side Wd 1 . The visually brightest stars seem to be concentrated 
in the small core defined above, but, given that neither the less 
evolved supergiants (their progenitors) or the WR stars (their de- 
scendants) share this concentration, this may be a random effect. 
In view of this, we conclude that the data so far strongly point 
to a single, homogeneous age (though some minimal spread, 
~ 1 Myr, is allowed) for the whole cluster. 

5. Conclusions 

Using new classification criteria, developed in Appendix A, we 
provide spectral types for almost 60 OB supergiants in Wd 1. 
Only a small fraction of them had previous accurate classifica- 
tions. 

The stars leaving the cluster main sequence form a well pop- 
ulated clump at spectral types O9-B0, displaying a range of lu- 
minosities, which likely reflect a difference in initial rotational 
velocities and, perhaps, formation over a timespan. Sixteen lumi- 
nous supergiants, with spectral types in the range B1-B4 range 
span the temperature range between this clump and the three 
blue hypergiants that bridge the gap with the A/F hypergiants. 

After a careful analysis of colour excesses, we find no strong 
reason to suspect that the extinction law towards Wd 1 deviates 
strongly from the standard, and we attribute the slightly diver- 
gent E{B - V) values to the difficulty in obtaining precise B pho- 
tometry. The measured values of E{V-I) and E(V-R) show that 
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Fig. 10. /?-band finder (~ 6' x4') for all the stars classified as OB supergiants in this paper Other stars in more advanced evolutionary 
stages are also shown. The moderately evolved OB supergiants are shown as rectangles. The more evolved stars are shown as circles: 
blue circles are blue hypergiants and LBVs, red circles are A/F hypergiants and M supergiants, while magenta circles are WR stars. 
The inset (~ 70" x 70") shows the finder of the crowded cluster core. Many stars in this area (and also in the halo) have not been 
observed. Note that WR stains N, X and T are outside the area covered by this finder. 



the extinction varies strongly across the face of the cluster, with 
average values implying Ay - 10.8. This value is fully consistent 
with the A/fj = 1 . 1 3 ± 0.03 derived by Brandner et al. ( 2008) and 
a stand ard law. The value A/f^ - 0.96 found by Crowth er et alJ 
(l2006ah . though subject to much higher uncertainty, is also com- 
patible. 

The dereddened magnitudes of the OB supergiants support a 
distance to Wd 1 similar to that found by Crowther et al. (2006a) 
from the infrared magnitudes of WR stars, namely DM - 13.5 
(d - 5.0 kpc). With this distance modulus, the intrinsic mag- 
nitudes of OB supergiants correspond closely to their spectral 
types according to typical calibrations. The short er DM - 13.0 
{d = 4.2 kpc) found by KotheS-& Doughertvl (l2007h is dis- 
favoured (though their error bars re ach d - 4.9 kpc). T he 
DM = 12.75 {d = 3.5 kpc) found bv iBrandner et al.1 (I2OO8I) is 
strongly disfavoured, as it would mean that all our la supergiants 
would have intrinsic magnitudes typical of lab supergiants, im- 
plying masses % 25 Mq, very difficult to reconcile with their age 
derivation of 3.6 Myr. We suspect this to be due to the sy stemat- 
ically younger ages that pre-MS isochrones indicate (cf. i Navlorl 
I2009h . 

The picture drawn from the distribution in spectral types 
seems compatible with a single burst of star formation and 



the predictions of evolutionary models until the stars start to 
move quickly towards lower temperatures. Further studies, in- 
cluding the use of the population of unevolved stars to nail 
down the cluster age and detailed analysis of the evolved popula- 
tion (which will provide accurate stellar parameters and element 
abundances), will undoubtedly result in very firm constraints on 
the evolutionary tracks followed by massive stars after they leave 
the main sequence. The high binary fr a ction found amongst mas- 
sive stars in Wd 1 dClark et allilOOl: iRitchie et al.ll2009ah will 
likely allow a better understanding of the role that binarity plays 
in this evolution. Furthermore, characterisation of binary prop- 
erties (frequency, mass ratio, period distribution, etc.) in the ho- 
mogeneous massive star population of Wd 1 can provide strong 
constraints on the formation mechanism of massive stars in such 
dense environments^, as the different for mation mechanisms 
propo sed will predict different outcomes (IZinnecker & York3 
l2007i for a review). Long-term spectroscopic monitoring of bi- 



^ Note, however, that at an age ~ 5 Myr, Wd 1 must have already un- 
dergone dynamical evolution. Indeed. ISrandner et alj j2008i) find some 
evidence for mass segregation, which could partially erase the signature 
of the formation mechanism. 
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naries in the cluster is underway (iRitchie et al.ll2009ah . and will 
deliver some of the observational characteristics. 
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Appendix A: Spectral classification of OB 
supergiants in the red 

A.1. Historical introduction 

OB stars have spectra characterised by strong B aimer and He 
lines. Metallic lines are few and weak, except in the range 09- 
B2, especially in main sequence stars. The blue range concen- 
trates most of the salient spectral features and spectral classifi- 
cation has b een developed using features in the /t/t 3900- 4700 
range (e.g.. i WalbornI Il97l IWalborn & Fitzpatricklll99ai) . The 
red part of the spectrum contains fewer features and has gen- 
erally been ignored, except for Ha, considered a good tracer of 
mass loss. An overview of features in the AA 5400 - 6600A range 
is presented by Walborn (1980), together with a spectral atlas 
based on photographic spectrograms. High-quality digital spec- 
trograms in this ra nge of a few selected objects are presented in 
IWalborn & Howar th (2000). 

Further to the red, atmospheric features become a compli- 
cation. There is, however, an important windo w, free from at- 
mosph eric features, in the AA 8300-8900 range. lAndrillat et"an 
(Il995h present spectra of early-type stars in the AA 8375 - 
8770 range, concluding that supergiants can be distinguished 
from main sequence stars from the shape of the Paschen lines. 
According to these authors, stars earlier than B3 cannot be clas- 
sified in this region. Higher resolution spectra were presented by 
iMunari & Tomasellal (Il999l) . Their spectral atlas covers the AA 
8480-8750 ran ge ai R = 2Q 000, but little attention is paid to 
early type stars. iMunari & Tomasellal (1 19991) find some degen- 
eracy between luminosity class and rotational speed and con- 
clude that O-type stars cannot be classified at all with this limited 
range. 

ICaron et"an (l2003l) studied quantitative classification crite- 
ria, using intermediate-resolution spectra in the AA 8400 - 8900 
range for stars with spectral types between 09 and B5. They 
found that the lines are always broader in dwarfs than in super- 
giants. In addition, the equivalent widths of most lines increase 
with advancing spectral type. Based on these findings, they de- 
veloped a method for approximate classification based on line 
ratios. This approach was followed in C05, where the grid was 
extended by using synthetic spectra and the validity of the classi- 
fications was checked against spectra of stars of known spectral 
types (see Appendix A in lC05h . 

Quantitative analysis of features in the red region have 
hardly been reported. The only exceptions are Cii 6578A, 
Si n 6347A and Ne i 6402A, included by Lennon et al.' ('1993) in 
their analysis of fea tures seen in B-type supergiants. In addition, 
iDavies et al.l ( l2005h present an atlas of luminous supergiants in 
this range. 



A.2. Data used 

In an attempt to improve our knowledge of the behaviour of fea- 
tures in the red and far red regions, we have collected spec- 
tra of OB stars i n this range from a n umber of sources. The 
STELIB library dLeBorgneet al.ll2003h contains stellar spec- 
tra at 7? ~ 2000. The data cover the whole optical range, but 
the number of OB su pergiants observed is small. The library by 
ICenarro et all (12001') covers the .1.18350-8900 range at better 
spectral resolution, but, being conceived to study the Ca ii triplet, 
contains few early-type stars. The more recent MILES library 
dSanchez-Blazquez et alj|2006h . which includes more early-type 
stars, only extends until 7500A. 

Th e main sources o f data have been the Indo-US li- 
brary dValdes et al.ll2004 h. which includes spectra in the range 
/1/13460 - 9460 at intermediate resolution and the UVES Paranal 
Observatory Project (ESO DDT Program ID 266.D-5655) spec- 
tra (Bagnulo et al. 2003). The UVES POP spectra have excellent 
resolution 7? 80 000 and SNR. Unfortunately, as the data were 
taken as part of a poor weather backup programme, the red sec- 
tions are generally affected by very strong atmospheric features. 
Moreover, the spectra present a gap in the /1/18540-8660 inter- 
val, due to the configurations used. 

In all, the spectra represent a rather complete, if heteroge- 
neous, sample of B-type luminous stars. The O subtypes are 
more sparsely sampled. 

A.3. Results 

In spite of the scarcity of spectral features in the red spectra of 
OB stars, several criteria can be used for spectral classification. 
Unfortunately, most of these criteria are sensitive to both spec- 
tral type and luminosity, but by considering the whole spectrum 
in the /1/16000 - 9000 range and combining several criteria, we 
can narrow the spectral range for a given supergiant spectrum 
generally down to half a spectral subtype. The classification of 
dwarf stars is much less accurate. In what follows, we will con- 
centrate on the features used to classify the supergiants seen in 
Wd 1. 

A.3.1. Classification in tine /-band 

Criteria for the class ification of 09-B 5 stars in this range have 
been laid down by ICaron etaTI (1200 3). They base their clas- 
sification in the relative strengths of Paschen lines. The best 
temperature criteria seem to be the EW(Pa 16+Ciii)/EW(Pa 11) 
or EW(Pal5)/EW(Pall) ratios, even though Hei 8850A 
is blended with Pa 1 1 at moderate resolutions. The ratio 
FWHM(Pa 13)/FWHM(Hei 8779 A) is a good luminosity crite- 
rion at the resolution they use. Unfortunately, in the spectra of 
stars in Wd 1, suffering from heavy extinction, there seems to 
be a complex system of diffuse interstellar bands. The 8620A 
band is sharp and narrow. Munari et al. (2008) has reported an 
excellent correlation between its strength and reddening up to 
rather high reddenings, at least E(B - V) ~ 2. However, at the 
much higher reddening of Wd 1, this feature is contaminated 
by neighbouring ones. There is a broad, relatively shallow fea- 
ture, perhaps similar in shape to the 4430A DIB, centred around 
8648 A. This feature is blended with a Hei 865 3 A line, w eakly 
seen in many spectra of B-type stars in lCaron et al. I (l200l . 

In addition, there appears to be a global depression of the 
continuum in the whole range 8540-8660A. The actual nature 
of this depression is difficult to decide, but it has a strong ef- 
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Fig.A.l. Luminosity sequence at B3. The unlabelled spectrum 
at the top is W70 (B3 la). All the other spectra have been taken 
from the Indo-US library. Note the strengthening and sharpening 
of He I and N i lines with increasing luminosity and the changing 
Pa I8/O1 8446 A ratio. The spectrum of W70 is affected by the 
strong interstellar features mentioned in the text. 
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Fig.A.2. Spectral seq uence of high-lum inosity supergiants. 
Spectra are taken from Cenarro et al. (1200 lb and have only mod- 
erate resolution and occasionally low SNR. Note the increase of 
N I lines with advancing spectral type and the sharp maximum 
of the He i around B2. Note also that the lines marked above the 
top spectrum are those seen in B-type spectra and do not match 
those on the AO la spectrum. 



feet on continuum normalisation. Treating it as a single feature 
produces an awkward shape for the continuum and results in a 
wrong sequence of strengths for the Paschen lines. In view of 
this, for the purposes of this paper, we have normalised the spec- 
trum as if the continuum around 861 OA were not affected by 
this broad depression (effectively assuming that there are two 
such depressions, one on each side of 861 OA) and we have nor- 
malised the spectra taking reference points in the continuum at 
8520, 8610 and 8690A. This procedure results in a decrease of 
the strength of the Paschen lines as we move up the sequence and 
it is hence taken to be approximately correct. However, it is only 
an approximation and we assume that the equivalent widths of 
Pa 13, 14 and 15 cannot be measured with accuracy. Because of 
this, we are forced to try to define a morphological classification 
scheme. 

Several Hei lines are seen in the /-band spectra of B-type 
stars. Though never near as prominent as the Paschen lines, 
they are much better defined in supergiants, showing maximum 
strength at B2S In dwarf stars. He i lines are always weak and 
only seen in the 08-B2 range (perhaps visible in high-SNR spec- 
tra of B3 V slow rotators). In giants, they are visible up to B4-5, 
depending on Vrot- In supergiants, they are seen down to B6 or 
even later if the SNR is very good. 

The O I 8446A line has a strong dependence on both spectral 
type and luminosity. As shown in Fig. lA.ll at a given spectral 

^ The Hei 858 2.7 & 8776.7 A lines were first described by 
Ivan Helder] il972h . who discusses their sensitivity to luminosity 



type, the ratio with the neighbouring Pa 18 (8437 A) decreases 
slowly with increasing luminosity. Figure [A2l on the other hand, 
shows how the same ratio increases with spectral type at a given 
luminosity. To make things worse, the shape of the line depends 
strongly on the rotational velocity. This can be easily seen by 
comparing the spectrum of HD 120315 (B3 V, Vrot - 150kms_^ 
and H P 160762 (B3IV and a very slow rotator; see lAbt et atl 
2002). In any case, O i 8446A is first seen at B2 in dwarf stars 
and starts to appear on the wing of Pa 18 at B 1.5 in supergiants. 

An important diagnostic for B3 is the appearance of Ni 
lines. They are weak at this spectral type, but quickly become 
stronger at later types, becoming prominent in supergiants after 
B5. These lines are not seen in dwarf stars. They are not seen 
in most spectra of luminosity class III objects, even the high- 
resolution UVES spectra. They are clearly visible in the few lu- 
minosity class II spectra available (including the B3 II star in 
Fig. lA.lb . and so it is reasonable to expect that they may appear 
weakly in some giants wit h low Vrot. In fact, they are weakly vis- 
ible in the spectrum from lAndrillat et al.l (1 19951) of H P 186568 
(B8III), which is a very slow rotator (vrot = 15kms~': lAbt et all 
2002). 

As discussed in |C05L O-type stars can be easily identified 
by the presence of Heii 8238A. Unfortunately, this line lies in 
the middle of a strong atmospheric band. When telluric features 
are accurately removed. Hen 8238A can be readily identified 
(Fig. IA.3b . Then it is visible up to at least 09.5, though it may 
be difficult to locate in fast rotating dwarfs. When the telluric 
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Fig. A.3. /-band spectra of O- and early B-type stars. These spec- 
tra from the Indo-US library have had the telluric features re- 
moved and allow the study of the behaviour of He ii 8238 A. Note 
how this line clearly identifies O-type stars. Note also the much 
higher number of Paschen lines seen in supergiants at a given 
spectral type. Finally, the C iii 8500A line is an excellent indica- 
tor of late O-type stars. 



features cannot be removed, as in the UVES POP spectra, it is 
almost impossible to see, though it may be guessed in the spectra 
of early and mid O-type stars. Based on the few spectra avail- 
able, this line may show some dependence on spectral type, be- 
ing weaker in 08-9 stars. No strong dependence on luminosity 
class can be identified. 

The second important marker of O-type stars is C iii 8500A. 
This line almost coincides with Pa 16 (8502 A) in wavelength. In 
06-7 dwarfs, where Pa 16 is hardly seen, it appears as a weak 
narrow line (Fig. IA.3I ). In O8-B0 stars it can be noticed as a 
important strengthening of Pa 16 in comparison with Pa 15 and 
Pa 17. This strengthening stops being noticeable around BO. 5, in 
line with the behaviour of C iii lines in the blue. As other C lines, 
this might be an unreliable tracer in stars with CNO anomalies, 
though we lack appropriate spectra to test this issue. 

A.3.2. Classification in tine red region 

Spectral features in the red region are scarce in OB stars. 
Many dwarfs show only Ha and the Hei 6678 & 7065A lines. 
Supergiants present a few other interesting features, notably in 
B-types. 

There are two He ii lines present in O-type spectra (a third 
one. He ii 6560A, cannot be separated from Ha). Unfortunately, 
He II 6527A coincides with a strong DIB and cannot be used for 
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Fig. A.4. Spectral features seen around Ha in O-type stars. The 
spectra are drawn from several sources and have quite different 
resolutions. The He 1 6678A line has not been marked to prevent 
confusion. 



classification. He ii 6683A starts to be seen as a weak inflection 
on the wing of He 1 6678 A in high-SNR spectra of 09 stars, and 
it becomes stronger at earlier spectral types. However, a clear de- 
pendence on spectral type cannot be drawn from the few spectra 
available (Fig.lA!4li . The ratio He 1 6678 A/He ii 6683 A cannot be 
used to determine spectral type. Indeed, from the four stars with 
spectral type 06-6.5, it would seem to depend more strongly on 
luminosity. 

A clear marker of high luminosity in O-type stars is the 
presence of the Si IV 6668, 6701 A lines in emission. These are 
"wind" or selective emission lines, due to non-LTE effects in 
an extended atmosphere. They are seen in luminous supergiants 
earlier than 09. We lack information on their behaviour in early 
O-types. The Siiv lines in the blue part of the spectrum go into 
emission at much earlier spectral types dWalborn et al.ll2002l) . 

Moving to later spectral types, we find another interesting 
"wind" feature, the C ii 723 1 , 7236 A doublet. Unfortunately, this 
feature lies in the middle of an atmospheric band, but, being in 
emission, is easily identifiable. This feature is visible between 
09.7 and B2. Its strength is clearly correlated to that of Ha emis- 



I. Negueruela et al.: OB supergiants in Westerlund 1 



15 




6550 6600 
Wavelength (A) 



Fig. A.5. Spectral features seen around Ha in B-type supergiants 
stars. The spectra are drawn from several sources and have quite 
different resolutions. Spectra from the UVES POP show strong 
atmospheric features. The correlation between the intensity of 
the C II 6578, 6582A doublet and spectral type is loose because 
of CNO anomalies, as proved by the strong anticorrelation with 
the strength of N ii lines at a given spectral type. 



sion, though there may be some dependence on spectral type, as 
well, with stronger features being seen in the BO. 5-1 range. It is 
restricted to stars with heavy mass loss, and so generally only 
appears in luminosity class la supergiants. 

The be haviour of Hg is similar. This line is a good tracer of 
mass loss (lLeithereilll9"88l) . It tends to be in emission in lumi- 
nosity class la supergiants (as shown in Fig. IA.5I ) and appears 
normally in absorption in luminosity class lb stars. Luminosity 
class lab supergiants occupy a middle ground, with Ha weakly 
in absorption or neutralised on many occasions, but in emission 
or strong absorption in other cases. There a re many exceptions 
to this general rule (see more examples in iDavies et al.l 120051 
Appendix A). In many case s, variable Ha emission is taken as 
evidence for binarity (e.g.. lThalleij|1997h . but monitoring of a 
few luminous OB supergiants has sho wn that im portant changes 
can happen on very short timescales (iMorel e t al. 2004). This 
behaviour has also been observed amongst the brightest OB su- 
pergiants in Wd 1 (Clark et al. 2010). 

The shape of Ha in emission has some dependence on spec- 
tral type. In B-type supergiants, it appears as a P-Cygni pro- 
file, sometimes superimposed on a broader plinth. The absorp- 
tion trough becomes weaker around B 1 -2 and tends to disappear 
completely. BO-0.5 supergiants generally show a strong single 
peak. O-type supergiants tend to have more complex features, 
with more than one peak. This is complicated by the presence of 
the He ii 6560A line in absorption. 

Emission is weaker for late-B supergiants, which gener- 
ally have weak P-Cygni profiles, but there are also many ex- 



ceptions. For example, the BSIa supergiant HD 199478 can 
show a strong Ha single peak emission, though this line is 
variable (Markova & Valchev 2000). The most luminous hy- 
pergiants may also display strong symmetric peaks. For exam- 
ple HD 32034 (B8-9Ia+ ) or H P 269781 (A2Iae) in the cata- 
logue of lLe BorgneetaLl (l2003l) do, but HD 33579 (A3 Ia+) and 
HD 223385 (A3 la^) display a weak absorption trough super- 
imposed on the strong emission peak (Note that the extremely 
high luminosity of all these stars, except HD 223385, can be 
confirmed, as they are members of the LMC). 

The Oi 7774A triplet is a very good luminosity indicator. 
Calibratio ns of its strength against l uminosity exist for A-G 
stars (e.g.. lAre llano Ferro et al. I l200l . but its behaviour has not 
been tabulated in B-type stars. The Unes are barely visible in B 1 
supergiants at moderate resolution, but quickly grow in strength 
with advancing spectral type. They are seen in dwarfs from B2- 
3, depending on rotational velocity. 

The C n 6578, 6582A doublet is generally not seen in MS 
B-type stars, where it is very weak and easily masked by high 
rotation or low SNR. It is present in B2-B5 stars of moderate lu- 
minosity and in most B-type supergiants (Walborn 1980). As it 
is clearly seen in many luminosity class IV stars, we suspect that 
it may be seen in MS stars with very low Vrot- At high SNR, it is 
weakly visible in BO. 5 supergiants. It clearl y appears in B0.7 su - 
pergiants and its strength peaks around B3 dLennon et alJll993h . 
The lines are visible until B8, and also in some B9 supergiants, 
with a clear dependance on luminosity. As seen in Fig. IA.5I 
their strength is not monotonically correlated to spectral type 
amongst luminosity class la supergiants. There seems to be some 
degree of anticorrelation between the strength of these features 
and those of N n, suggesting that differences in CNO abundances 
play a role in their relative strength, as is seen in features further 
to the blue d Walborn & HowartljhoOOl) . 

There are two Nii lines in the range showed in Fig. IA.5I 
which are only visible in supergiants. N ii 6482A appears weakly 
in Bl supergiants, is rather strong at B2 and becomes weak 
again at B3, disappearing at later spectral types. Its dependence 
on temperature is very similar to that of Nn lines in the blue 
range, such as Nn 4601 A, tabulated bv iLennon et al. I ([T99I . 
Nn 661 OA is rather we aker, and has t he same temperature de- 
pendence. As shown by Lennon ( 1993), these two lines are only 
prominent in B2 la supergiants. However, N-enhanced stars of 
similar spectral types may also show relatively strong features 
(Fig. IA.5b . Another Nn line behaving similarly is Nii 6380A 
(Walborn 1980). 

The Ne i 6402A line i s seen in seve r al of t he spectra. Its be- 
haviour was tabulated bv iLennon et al. l (fT99l . The line is visi- 
ble from Bl to B9, and strongest at B4-5, but it does not seen to 
depend on luminosity subclass amongst supergiants. A weaker 
line, Ne i 6506A, is seen in the spectrum of HR Car, which is 
classified B 2 dNota et al.lll997l). The Sin 6347A line was also 
tabulated by Lenn on et al .I (ll993h . It is hardly visible at B2, but 
increases its strength very quickly, becoming prominent at B3 
and growing in strength towards later spectral types. The neigh- 
bouring Sin 6371 A behaves s imilarly. Features further to the 
blue are discussed bv lWalbornl(IT980.) . 



A.4. Summary and outlook 

The results presented here show that classification of early-type 
stars is possible using only red spectra. A summary of features 
useful for classification is presented in Table A.l. The classi- 
fications are not as accurate as those using blue spectra, but 
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likely comparable to classifications based on K-hand spectra. 
Combination of red spectra with near-lR spectra offers good 
prospects to classify rather accurately obscured early-type stars. 
The limited spectral range offered by GAIA does not allow ac- 
curate spectral c lassification, but the s i tuation is not as desperate 
as suggested bylM unari & Tomasellal ([T999h . As seen in Fig. |6] 
even O-type stars can be approximately classified in that range, 
at high resolution and signal-to-noise ratio. 

Further improvement on classification and a more thorough 
understanding of some of the issues raised here may be obtained 
by a systematic study of MK standards and stars with accurate 
spectral types. 
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I. Negueruela et al.: OB supergiants in Westerlund 1 
Table 1. Sample of OB supergiants in Wd 1, with BVRI photometry from lCOSi (when available) and spectral types'". 
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ID 


RA (J2000) 


Dec (J2000) 


B 


V 


R 


/ 


Spectral 
Type 


W2a 


16 46 59.71 


-45 50 51.1 


20.4 


16.69 


14.23 


11.73 


B2Ia 


W5 


16 47 02.97 


-45 50 19.5 


21.4 


17.49 


14.98 


12.48 


B0.5Ia+ 


W6a 


16 47 03.04 


-45 50 23.6 


22.2 


18.41 


15.80 


13.16 


B0.5Iab 


W6b" 


16 47 02.93 


-45 50 22.3 


23.6 


20.20 


17.91 


15.25 


09.5 III 


W7 


16 47 03.62 


-45 50 14.2 


20.0 


15.57 


12.73 


9.99 


B5 Ia(+' 


W8b 


16 47 04.95 


-45 50 26.7 


- 


- 


- 


- 


B1.5Ia 


WIO 


16 47 03.32 


-45 50 34.7 


- 


- 


- 


- 


B0.5 1+OB 


Wll 


16 47 02.23 


-45 50 47.0 


21.2 


17.15 


14.52 


11.91 


B2Ia 


W13 


16 47 06.45 


-45 50 26.0 


21.1 


17.19 


14.63 


12.06 


B0.5Ia++OB 


W15 


16 47 06.63 


-45 50 29.7 


22.8 


18.96 


16.38 


13.75 


09 lb 


W17" 


16 47 06.25 


-45 50 49.2 


22.7 


18.87 


16.19 


13.56 


09 lab 


W18 


16 47 05.71 


-45 50 50.5 


21.2 


17.32 


14.81 


12.27 


B0.5 la 


W19 


16 47 04.86 


-45 50 59.1 


22.6 


18.22 


15.21 


12.37 


Blla 


W21 


16 47 01.10 


-45 51 13.6 


22.5 


18.41 


15.56 


12.74 


B0.5Ia 


W23a 


16 47 02.57 


-45 51 08.7 


22.1 


17.85 


14.91 


12.07 


B2Ia+B I? 


W24 


16 47 02.15 


-45 51 12.4 


23.0 


18.71 


15.96 


13.24 


09 lab 


W25* 


16 47 05.78 


-45 50 33.3 


21.9 


17.85 


15.22 


12.61 


09Iab 


W27* 


16 47 05.15 


-45 50 41.3 


21.5 


17.94 


15.35 


12.80 


09I+OI 


W28 


16 47 04.66 


-45 50 38.4 


20.9 


16.87 


14.26 


11.64 


B2Ia 


W29 


16 47 04.41 


-45 50 39.8 


22.6 


18.66 


16.02 


13.38 


09Ib 


W30 


16 47 04.11 


-45 50 39.0 


22.4 


18.45 


15.80 


13.20 


O+O 


W31 


16 47 03.78 


-45 50 40.4 


- 


- 


- 


- 


BOI+OB 


W33 


16 47 04.12 


-45 50 48.3 


20.0 


15.61 


12.78 


10.04 


B5 Ia'+' 


W34 


16 47 04.39 


-45 50 47.2 


22.1 


18.15 


15.40 


12.69 


BO la 


W35 


16 47 04.20 


-45 50 53.5 


22.7 


18.59 


16.00 


13.31 


09Iab 


W37" 


16 47 06.01 


-45 50 47.4 


22.8 


19.11 


16.40 


13.65 


09Ib 


W38" 


16 47 02.86 


-45 50 46.0 


23.2 


19.10 


16.47 


13.81 


09 lab 


W41 


16 47 02.70 


-45 50 56.9 


21.3 


17.87 


15.39 


12.78 


09 lab 


W42a 


16 47 03.25 


-45 50 52.1 


- 


- 


- 


- 


B9 Ia'+> 


W43a 


16 47 03.54 


-45 50 57.3 


22.8 


18.05 


15.22 


12.26 


BO la 


W43b 


16 47 03.52 


-45 50 56.5 


- 


- 


- 


- 


Blla 


W43c 


16 47 03.76 


-45 50 58.3 


20.4 


18.35 


16.18 


13.66 


09Ib 


W46a 


16 47 03.91 


-45 51 19.5 


23.0 


18.55 


15.46 


12.46 


Blla 


W46b 


16 47 03.61 


-45 51 20.0 


- 


- 


- 


- 


09.5 lb 


W49" 


16 47 01.90 


-45 50 31.5 


22.6 


18.76 


16.30 


13.80 


BO lab 


W50b** 


16 47 01.17 


-45 50 26.7 


22.8 


19.66 


17.21 


14.69 


09 III 


W52 


16 47 01.84 


-45 51 29.2 


21.8 


17.48 


14.68 


11.94 


B 1.5 la 


W54 


16 47 03.06 


-45 51 30.5 


- 


- 


- 


- 


B0.5Iab 


W55" 


16 46 58.40 


-45 51 31.2 


21.6 


17.67 


15.25 


12.67 


BO la 


W56b** 


16 46 58.85 


-45 51 45.8 


22.8 


18.88 


16.36 


13.76 


09.5 lb 


W57a 


16 47 01.35 


-45 51 45.6 


20.7 


16.54 


13.83 


11.13 


B4Ia 


W60 


16 47 04.13 


-45 51 52.1 


22.8 


18.50 


15.96 


13.28 


BO lab 


W61a 


16 47 02.29 


-45 51 41.6 


21.2 


17.16 


14.62 


12.01 


B0.5Ia 


W61b 


16 47 02.56 


-45 51 41.6 


22.7 


18.59 


16.00 


13.31 


09.5 lab 


W62a 


16 47 02.51 


-45 51 37.9 


- 


- 


- 


- 


BO. 5 lb 


W63a 


16 47 03.39 


-45 51 57.7 


22.6 


18.56 


16.20 


13.68 


BO lab 


W65" 


16 47 03.89 


-45 51 46.3 


22.9 


18.73 


16.27 


13.68 


09Ib 


W70 


16 47 09.36 


-45 50 49.6 


21.2 


16.88 


14.10 


11.29 


B3Ia 


W71 


16 47 08.44 


-45 50 49.3 


21.5 


17.01 


14.06 


11.16 


B2.5 la 


W74 


16 47 07.08 


-45 50 13.1 










09.5 lab 


W78 


16 47 01.54 


-45 49 57.8 


21.0 


17.06 


14.54 


12.04 


Blla 


W84" 


16 46 59.03 


-45 50 28.2 


21.3 


17.82 


15.60 


13.63 


09.5 lb 


W86" 


16 46 57.15 


-45 50 09.9 


22.9 


18.76 


16.43 


14.00 


09.5 lb 


W228b 


16 46 58.05 


-45 53 01.0 










09Ib 


W232 


16 47 01.41 


-45 52 34.9 


21.3 


17.53 


15.25 


12.85 


BO lab 


W238 


16 47 04.41 


-45 52 27.6 


21.4 


17.47 


14.98 


12.45 


Bllab 


W373" 


16 46 57.71 


-45 53 20.1 










BO lab 



" Objects marked with "*" have classifications based on EMMI low-resolution (grism #6) spectra. Their spectral types are therefore less 
accurate. 

Objects marked with "**" have classifications based only on FLAMES spectra covering 8480-8920A. Their spectral types are tentative. 
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Table A.l. Summary of features useful for spectral classification between 04 and B9 in the 6300-9000A range. 







to 


11 Id-Aiiii mil 




Paschen lines 


04 


B9 


>B9 


Number of lines seen depends on temperature and luminosity 


He I 6678, 7065A 








Shape of lines correlates with luminosity 


04 


~B7V, B91 


B2 




He I 8847, 8777, 8583A 


08-9 


B3V, B61 


B2 


strength depends on luminosity 


Hen6683A 


04 


07V, 091 


04 


blended with He i 6678A 


He II 8283 A 


04 


091-V 


04 


inside telluric band 


Cm 8500A 


-06 


B0.5 


09 


lower limit not well determined 


Cn6578, 6582A 


B0.7 


B6111, B81 


B3 


rarely seen in dwarfs 


Cn 723 1,7236 A 


09.7I 


B21 


Bl 


emission features only in luminous SGs 










inside telluric band 


Oi8446A 


B2VB1.51 


B9 


>B9 


strong dependence on luminosity 


1 7774A 


B3VB1I 


B9 


>B9 


strong dependence on luminosity 


Si II 6347, 6371 A 


B21 


B9 


B8 


tabulated bv Lennon et al. ("1993) in SGs 


Si IV 6668, 6701 A 




09 


? 


emission lines in SGs; lower limit not determined 


N I 8680-86, 8703-18A 


B3 1-11 


B9 1-11 


>B9 


not seen in 111-V 


Nil 6482, 66IOA 


Bll 


B31 


B2 


only strong in la SGs 


Ne I 6402 A 


Bll 


B91 


B4-5 


tabulated bv Lennon et al. ("1993^ in SGs 



